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Abstract. The bulk viscosity of cold, dense three-flavor quark matter is studied as 
a function of temperature and the amplitude of density oscillations. The study is 
also extended to the case of two different types of anharmonic oscillations of density. 
We point several qualitative effects due to the anharmonicity, although quantitatively 
^ , they appear to be relatively small. We also find that, in most regions of the parameter 

pH ' space, with the exception of the case of a very large amplitude of density oscillations 

(i.e. 10% and above), nonlinear effects and anharmonicity have a small effect on the 
, interplay of the nonleptonic and semileptonic processes in the bulk viscosity. 

in 

Q ■ 1. Introduction 

(N 

Cold, dense quark matter is one of possible states of baryonic matter formed at very 
high densities. While there is little doubt that such quark matter can be formed in 
principle, the value of the critical density needed remains unknown. This is one of the 
reasons why it is not settled whether such form of matter can exist inside neutron stars. 
' Indeed, the interior region of neutron stars is the most likely place to find dense quark 
matter. The corresponding densities reach up to about 10 times the nuclear saturation 
density, while the temperatures remain moderately low, of the order of 1 MeV or less. 

A recent report on a precision measurement of the mass of the binary millisecond 
pulsar J1614-2230 [1] seems to strongly constraint the possibility of stellar dense quark 
matter [2]. Nevertheless, it is hard to rule out quark matter solely based on the 
measurement of such a global property of a star such as its mass.. More informative 
probes of the state of matter at the highest densities are likely to be based on stellar 
characteristics determined by the transport properties of matter, various emission rates, 
as well as certain thermodynamic properties, which are sensitive to the spectrum of 
quasiparticles at the Fermi surface and which, therefore, can provide a deeper insight 
into the microscopic nature of dense matter. 

In this paper, we study one of such transport characteristics of dense quark matter, 
the bulk viscosity. In general, viscosity is one of possible mechanisms responsible for 
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damping of the so-called r-mode instabilities in compact stars [21 IH El El E] • The emission 
of gravitational waves tends to drive a differential collective motion of the stellar fluid 
in the form of r-modes |H1 El EDI [H] • Such modes in turn increase the stellar quadrupole 
moment and feed back to result in even stronger gravitational emission. If not damped 
by dissipative processes, the resonant growth of the r-modes can eventually lead to a 
breakup of the star. 

The bulk viscosity in the normal phase of three-flavor quark matter is usually 
dominated by nonleptonic weak processes [31 [71 [121 ESI [HI [ISl [iSl [13 [H], shown in 
Figs. [U (a) and (b). Under certain conditions, semileptonic processes, see Figs. [1] (c)- 
(f), can lead to an order of magnitude increase of the viscosity as a result of their subtle 
interplay with the nonleptonic processes [191 [20] • 

When the temperature is sufficiently low or when the magnitude of the density 
oscillation is sufficiently large, the rates of the weak processes may get substantial 
nonlinear dependence on the parameter 5/ij/T, where are the chemical potential 
imbalances that measure the departure of quark matter from /3 equilibrium. In the 
context of the nonleptonic processes in strange quark matter, such effects have already 
been studied [31 [151 IlH] (see also Ref. [21]). Here we extend the analysis to include also 
the semileptonic processes with the nonlinear corrections in the rates. 

When the nonlinear regime is realized, it may be reasonable to expect 
additional complications in the dynamics responsible for the energy dissipation of the 
hydrodynamic flow. In this paper we address one of such complications that is associated 
with the anharmonicity of the density oscillations driven by the r-modes. To the best of 
our knowledge, previously this has not been addressed in the literature. In studies of the 
bulk viscosity of stellar matter, one commonly assumes that the oscillations are perfectly 
harmonic. Of course, this is justifled in the linear regime, when different harmonics do 
not interfere. 

In general, large magnitude density oscillations can hardly remain perfectly 
harmonic because higher harmonics can be generated by the nonlinearity. As discussed 
in Refs. [221 1231 1211 [2S1 [23 [2Z] , the nonlinear regime in compact stars may be responsible 
for several qualitatively new features in the dynamics, e.g., coupling of eigenmodes, a 
non-sinusoidal shape of oscillations, various resonance phenomena, and the formation 
of shocks and turbulence. At the fundamental level, the nonlinear effects are the 
consequence of the equation of state of matter and the Einstein equations, e.g., see 




Figure 1. Diagrammatic representation of the weak processes that contribute to the 
bulk viscosity of dense quark matter. 



Bulk viscosity in the nonlinear and anharmonic regime 



3 



Ref . • Usually, they play a profound role when the interaction energy is comparable 
to the energies of modes. However, the nonlinearity is present in principle even at 
arbitrarily small amplitudes of oscillations. Then, one of the natural questions about 
the nonlinear regime is: How does the associated non-sinusoidal shape of oscillations 
(i.e. the anharmonicity due to mode coupling) change the bulk viscosity? It is the goal 
of this paper to make a first attempt in exploring this issue. Our study will be limited 
to the case when the density fluctuations are much smaller than the equilibrium density, 
although they can be much larger than the temperature. In view of this limitation, we 
cannot address many features of dynamics that may develop in the extreme nonlinear 
regime. Also, instead of considering realistic nonlinear oscillations found in stellar 
simulations [221 ISSl 1211 [251 |26l [27], we will simply model the shape of density oscillations 
by using classical solutions to an anharmonic oscillator problem. Such a toy-model study 
will hopefully lead to a better understanding of the problem and will trigger further 
developments in the future. 

The rest of the paper is organized as follows. In the next section, we present 
the general formalism for calculating the bulk viscosity in dense quark matter with 
multiple active weak processes and anharmonicity included. In Sec. [3l we start our 
study of the bulk viscosity in the harmonic regime of density oscillations. Our results 
qualitatively reproduce the earlier results of Refs. [3l [18] and extend them to include 
additional effects due to semileptonic weak processes. In the same section, we also 
comment on the linear regime when the expression for the bulk viscosity can be obtained 
analytically. The anharmonic regime is studied in Sec. [4|for two different types of density 
oscillations, modeled by solutions to a classical anharmonic oscillator with cubic and 
quartic potentials, respectively. A brief discussion of our result is given in Sec. [5] General 
periodic solutions in the case of classical anharmonic oscillators with two different types 



of potentials are presented in analytical form in Appendix A 
2. Formalism 

One can calculate the bulk viscosity C, under conditions realized in stars by comparing 
the hydrodynamic relation for the energy dissipation, averaged over one period r. 



(^diss) = -- / (V ■ vf dt ^ / (Snf dt. (1) 



^o'^ Jo 

with the thermodynamic relation for the mechanical work, counteracting the 
hydrodynamic flow, 

(^diss) = - / PVdt ~ / P6hdt. (2) 

r Jo riQT Jq 

The latter is given in terms of the instantaneous pressure P and the specific volume 
V = 1/n. In the above expressions, uq is the equilibrium density and 6n = n — uq is 
the density deviation from the equilibrium value. (Even when the magnitude of density 
oscillations is not vanishingly small, we will assume that \6n\ <^ no.) 
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If the magnitude of the density oscillations dn^ is vanishingly small, one may 
simulate the collective motion as a harmonic oscillation, 5n{t) = 6nQCOs{ujt). However, 
in a resonance regime, when large density oscillations develop, nonlinear effects may 
start to play an important role. In this paper, we study this possibility by simulating 
two different types of anharmonic density oscillations. The two types correspond to 
oscillators with cubic and quartic terms in the potential energy, which have different 
symmetry properties under 6n — —6n. The corresponding equations of motion read: 

5h + UqSu {1 + a5n) = 0, (Type I), (3) 

5n + uj^dn (1 + ^6n^) = 0, (Type II). (4) 

Note that the coupling constants a and (3 have the dimensions of an inverse density 
and an inverse density squared, respectively. It is convenient, therefore, to introduce 
the dimensionless parameters a* = aSno and /3* = /3(5no)^, which are given in terms 
of the amplitude of density oscillations Suq. [For asymmetric oscillations, described by 
Eq. ([3]), we assume that the amplitude is the maximum deviation from the equilibrium 
point.] Note that the parameters a* and /3* can be either positive or negative. General 
periodic solutions for both types of anharmonic oscillators can be given in terms of the 
Jacobi elliptic functions. The corresponding solutions are presented in Appendix A By 
substituting these exact solutions into Eq. ([1]) and making use of the result in Eq. ( ]A.4p . 
we derive 

(^.ss) = -^f^V^, (5) 



C = ^77r^Z / P^'^dt. (6) 



2 V ^0 

where the constant J-" for each type of solution is determined in terms of a* and /?*, see 
Eqs. (lA.Sp and flA.9p . respectively. As is easy to check, J-" — > 1 in the harmonic limit 
a* ^0 (Type I) or /3* ^ (Type II). 

By comparing Eqs. (j2]) and ([5]), we obtain the following expression for the bulk 
viscosity: 

2no 1 

When there is a departure from /3 equilibrium, the pressure can be given in terms of the 
instantaneous composition, 

dP 

P = P + —Sn + n{Ci-C2)6Xe + nCiSXs, (7) 
on 

where Xf, = Ue/n and X^ = Ug/n is the electron and strangeness fractions, and P is the 
pressure in equilibrium. The susceptibility functions Ci and C2 were defined in Ref . [19] . 
By taking into account that 5n is a periodic function, one finds that only the last two 
terms in the pressure ((Tj) contribute to the bulk viscosity ([6]), 

2t7^ 1 r 

C = 21, \2-r - / [(C*! - C2)SX, + C,6X,] 6n dt. (8) 

The instantaneous composition of quark matter is determined by the weak processes, 
shown in Fig. [H Taking all of them into account, we derive the following set of nonlinear 
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differential equations for the electron and strangeness fractions: 

(r. - rj + (r^ - r^) = x,5^^, g (^) 

+ A3 (5/i2 - f 



n ; — 

dt 



(9) 



n — — — = [Ta - Tb) + [Tc - Td) 



dt 



hs,,±r, (^py - x-M.j:xi {^f (10) 



j=0 ^ ' j=0 

where = /i^ — /i^, 6fi2 = fJ's — fi'u — /^e, and the notation for the A-rates are the same 
as in Refs. [I9l[20], 

64 

Ai = — G| sin^ 9, cos^ O^fxlT^ , (11) 

X2 = ^Glsm^e,f^smlT\ (12) 
4U7r 

17 

A3 = cos^ 9casiJ,dlJ'elJ'uT'^- (13) 

In Eqs. ([9]) and ffTOj) . all higher order corrections in powers of 6fii/T were taken into 
account, while higher order corrections in powers of 6fii/fi were neglected. This is the 
same approximation that was used in Ref. [18]. By calculating the semileptonic rates 
using the approach of Ref. [28], it is easy to check that the coefficients of the nonlinear 
terms are the same as in the nucleon direct Urea process [29], [30]: Xo = 1, Xi = t|^5 
and X2 = i7~4- The corresponding coefficients in the nonleptonic rates are Tq = 1 and 

^1 = ^ mmm- 

The functions that describe the deviation from equilibrium, Sfii, can be equivalently 
rewritten in terms of the electron and strangeness fractions: = Ci^+Bi6Xe+AiSXs, 
where coefficient functions Ai, Bi and Ci were defined in Ref. [19]. By making use of 
these relations and Eqs. ([9]) and ffTOl) . we derive the following self- consistent set of 
equations for the dimensionless quantities ui = 6fii/T (for i = 1,2), 

1 2 



i=o i=o 
2 

+ X,B,Y.X,{v2-yif'^'. (14) 
i=o 

In this study, 5n is a periodic function that describes anharmonic oscillations of either 
Type I or Type II, see Eqs. ^ and (jl]). We make use of the analytical results in 



Appendix A and solve Eq. ( [T4|) numerically. When the solutions for z/j (with i = 1,2) 



are available, one can invert the relations for in terms of 6n, 6Xe and 6Xs in order 
to determine the deviation of the electron and strangeness fractions, 

6Xe = Ge— + H^Pi + Jel^2, (15) 

no 
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5X, = G— + H,u^ + J,U2. (16) 
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where 



_ A1C2 — A2C1 _ B2C1 — B1C2 . „s 

~ A2B,-A,B2' ' ~ A2B,-A,B2' ^ ' 

TT _ TA2 _ TB2 , . 

A2B,-A,B2' ' A2B,~A,B2' ^ ' 

" ~A2B,-A,B2' ~ A2B,-A,B2 ^^^^ 

Finally, by making use of these results in Eq. ([8]), we can calculate the bulk viscosity, 

_ 2Tnl A2{C,-C2)-B2C, 1 r 
^ - uiiSnoyT A2B, - A,B2 r J, 

uii6no)^T A2Bi-A^B2 r ^ ^ ^ 

3. Harmonic oscillations 

In the limiting case of harmonic oscillations, the density deviations are described by 
6n = 6nQsm{uot) and the constant J-" in Eq. ([5]) is equal to 1. The numerical results 
for the bulk viscosity as a function of Suo/uq for several fixed values of temperature are 
shown in Fig. [2l The linear regime corresponds to small values of Srio/uQ , where the 
bulk viscosity saturates. It is also interesting to present the temperature dependence 
of the bulk viscosity. The corresponding plots for several fixed values of Suq/uq are 
shown in Fig. [3l As we see, with decreasing the temperature, the bulk viscosity 
eventually levels off. This is the outcome of reaching the nonlinear regime, and it 
would be absent if the linear approximation were used instead. In both Fig. [2] and 
Fig. [3], the results for two representative values of the period of oscillations, r = 0.1 s 
and r = 10~^ s, are shown. (In all calculations, we used the model parameters from 
Ref. [19] at n = 5pQ, where po — 0.15 fm~^ is the nuclear saturation density.) In most 
regions of the parameter space, our results qualitatively agree with earlier findings in 
Refs. [31 [T2I [131 EH [13 [m im [H] • The only notable difference occurs in Fig. [3] around 
the "semileptonic" hump (T ~ 1 MeV). This comes from the interplay of semileptonic 
weak processes with the more dominant nonleptonic ones [191 120] • In most of the studies, 
this is neglected because of a smallness of the semileptonic rates. 

Here it is appropriate to mention that, in application to stellar quark matter, the 
bulk viscosity may not be the only, or even the dominant mechanism responsible for 
damping of the r-mode instabilities. For example, at sufficiently low temperatures, the 
corresponding dissipative dynamics is known to be dominated by the shear viscosity 
[3 [32] . (For several representative studies of the shear viscosity in dense quark matter 
see, for example, Refs. [331 EH [35].) In this connection, our low-temperature results 
in Figs. [2] and [3] should be used only for the purpose of determining where exactly the 
damping by the shear viscosity takes over. 
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Figure 2. Bulk viscosity as a function of Sno/no for several fixed values of 
temperature, i.e. T = 10"'^ MeV (black solid line), T — 10~^ MeV (red long-dashed 
line), T = 10"^ MeV (blue dashed line), T = 10"^ MeV (green short-dashed line), 
T = 10"^ MeV (black dash-dotted line) and T = 1 MeV (red dash-dotted hne line). 
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Figure 3. Bulk viscosity as a function of temperature for several fixed values of the 
amplitude of density oscillations, i.e. 6no/nQ = 10"'^ (black solid line), Suq/uq = 10~^ 
(red long-dashed line), 6no/no = 10^^ (blue dashed line), 6no/no ~ 10^^ (green short- 
dashed line), and Sno/no = 10^^ (black dash-dotted line). 



In order to better understand the role of nonlinear terms in the weak rates, see 
Eqs. ([9]) and ( |T0|) . as well as their effect on the bulk viscosity, it is instructive to study 
the linear approximation. In this case, the expression for the viscosity can be derived 
analytically. As we shall see, this will be also helpful to elucidate the role of induced 
oscillations of Sfi2 = /^s — /^n — /^e when the semileptonic processes are formally switched 
off. 



3.1. Harmonic oscillations: linear approximation 

In the linear regime, Eqs. (|9]) and ( ITOl) for dimensionless functions Ui simplify down to 
di COS^? - fiUi - f2U2 + f3{t^2 - i^i), (21) 

(22) 



dui 
du' 



—— = 6/2 cos 1? - hiUi - h2V2 + h'i{v2 - ^l) 

dv 
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where = 27rt/r is the dimensionless time variable, di = CiSuo/ (Tuo) is the magnitude 
of the "driving force" , and the other coefficient functions are 

/i = Ai^, h, = \i^, (23) 

J2 = M , n-2 = A2 , (24) 

h = As-^, h, = Aa-^. (25) 

The general solution to this set of equations in the steady state regime is given by 

ui = xi cos + yi sin 1!}, (26) 
1^2 = X2 cos "d + y2 sin 'd, (27) 
where the coefficients satisfy the following set of algebraic equations: 

{h + h)xi + {f2-h)x2 + yi =d,, (28) 

-xi + ifi + h)yi + if2-h)y2 =0, (29) 
{hi + h3)xi + (/l2 - h3)x2 + y2 = d2, (30) 
-X2 + {hi + h3)yi + {h2 - h3)y2 = 0. (31) 

It is straightforward, although tedious to solve this set of equations. When the solution 
is available, the result for the bulk viscosity will follow from Eq. f l2U]) . i.e. 

^ _ Tnl A2{Ci - C2) - B2Ci ^^ ^ Tnl BiCi - Ai{Ci - €2) ^^ ^^^^ 
uSuq A2B1 — A1B2 uSuq A2B1 — A1B2 

It can be shown that this expression (with the appropriate solutions for xi and X2) 
coincides exactly with the result for the bulk viscosity, obtained in Ref. [T9] . 

3.2. Harmonic oscillations: nonleptonic contribution in linear approximation 

If one ignores the semileptonic processes (i.e. if one formally takes A2 = A3 = 0), the 
linearized equations f l2T|) and fl22|) for z/j's take the following form: 

^ = dicos?? - /ii/i, (33) 
dv 

—^ = d2COS-d-hiiyi. (34) 
d-d 

The explicit solution to this set of equations in the steady state regime is given by 

di 

i + n 

dihi 

and the corresponding expression for the bulk viscosity reads 

Tnl A2{Ci-C2)-B2Ci dji 



vi = ^ ^ (/i cosi? + sin'd) , (35) 
1 

U2 = d2 sin -d + ^ | (cos 'd — fi sin -d) , (36) 



Cn 



uduQ A2B1 -A1B2 1 + 

Tnl BiCi-Ai{Ci-C2) dihi ^ XiCf 

oo6no A2B1-A1B2 1 + + {XiAi/no)^' ^ ^ 
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where we used the relation fi/hi = A1/A2 to arrive at the final result. As expected, 
this agrees with the known result [3l [131 [HI [19] - 

It is interesting to notice that the final result for the bulk viscosity receives a nonzero 
contribution due to the oscillation of 1^2 = Since Sfj,2 = /Xg — /^u — /ig controls 

the imbalance of the rates in the semileptonic processes, shown in Figs. [1] (c) and [d), 
which are formally switched off in the approximation at hand, one might wonder why 
there should be such a contribution at all. The answer is quite simple. In absence of the 
semileptonic processes the electron fraction in quark matter cannot change. However, 
when the nonleptonic processes drive the oscillations of the strangeness composition, 
they inevitably induce the oscillations of 6fi2 = fJ's—fJ'u — fJ'e- Then, the latter contributes 
to the instantaneous pressure and, in turn, to the bulk viscosity. Interestingly, such a 
contribution due to the induced oscillation of 1^2 = 5 1x2 1 T were ignored in all previous 
studies [3l [121 [131 El IlSl [lEl [TTl [H]. Fortunately, the corresponding correction is 
quantitatively small. The reason for its smallness seems to be rooted in the "accidental" 
fact that one of the susceptibility functions, i?2, is inversely proportional to the square 
of the chemical potential of electrons (rather than quarks) and, thus, is considerably 
larger than all others p!9] . 

4. Bulk viscosity in anharmonic regime 

In this section we study the effect that anharmonic oscillations of the density have on 
the bulk viscosity of dense quark matter. 

4-1. Anharmonic oscillations of Type I 

Let us start by modeling the density oscillations of quark matter 5n{t) by a time 
dependent anharmonic function of Type I, which is a solution to Eq. Q with a 
fixed anharmonicity parameter a*. Before we proceed to the numerical results, it is 
important to notice that the corresponding oscillations are asymmetric with respect to 
the equilibrium point 5ncq = 0. For a* < 0, the density oscillations are larger in the 
direction of positive (5n, while for a* > 0, the oscillations are larger in the direction of 
negative 6n, see also Fig. lAli The cases of the positive and negative anharmonicity 
parameters are physically equivalent, however. Indeed, they are related by the following 
sign reversal symmetry: a — )■ —a and Sn — > —6n. Therefore, it is sufficient to study 
only one of them. For technical reasons, we choose a* < 0. 

Typical results for the bulk viscosity as a function of anharmonicity parameter a* 
are shown in the left panel of Fig. [Hfor the whole range of negative a*, i.e. — 1 < a* < 0, 
for which physically meaningful periodic solutions exist. We used the following values 
of the period and the amplitude of density oscillations: r = 10~^ s and 6nQ/no = 10~^, 
and plotted the results for several representative values of temperature. [It should 
be emphasized that the period r is related to the "bare" frequency Uq by a modified 
relation, see Eq. (1A.3I) .] In general, we find that the bulk viscosity decreases with 
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Figure 4. Bulk viscosity as a function of anharmonicity parameter of Type I 
(left panel) and Type II (right panel) for r = 10""^ s, dno/no = 10~^ and several 
representative values of temperature. 



increasing the degree of anharmonicity. This quahtative behavior may be understood 
as the result of an effective increase of the frequency of oscillations due to an admixture 
of higher harmonics. Quantitatively, however, the effect is rather small. Only a very 
large anharmonicity {a* ^ — 1) leads to a substantial decrease of the bulk viscosity. 

4-2. Anharmonic oscillations of Type II 

Anharmonic density oscillations of Type II are modeled by a function Sn{t), which 
is a solution to Eq. (jl]) with a fixed anharmonicity parameter Conceptually, this 
is a simpler case because the oscillations are symmetric about the equilibrium point 
6neq = 0. Unlike the case of Type I oscillations, there is no reversal symmetry here. As 
in the previous case, however, periodic solution exist only for a range of values of the 
the anharmonicity parameter, /3* > — 1. 

Numerical results for the bulk viscosity as a function of anharmonicity parameter 
/3* are shown in the right panel of Fig. HJ The qualitative dependence of the viscosity 
on the parameter /3* is somewhat different. While it decreases at large values of 
parameter there is a range of small negative values of where it slightly grows 
with increasing anharmonicity. Moreover, this feature seems to be rather general and 
especially pronounced in the nonlinear regime (small temperature). Just like in the case 
of Type I oscillations, the effects appear to be rather small. 

5. Discussion 

In this paper we studied the bulk viscosity of dense quark matter by taking into 
account the nonlinear dependence of the nonleptonic and semileptonic weak rates on 
the parameter 6fii/T, where are the chemical potential imbalances that control 
the departure of strange quark matter from (3 equilibrium. We reproduce the earlier 
obtained interplay of the nonleptonic and semileptonic processes, leading to an increase 
( "hump" ) of the viscosity in a narrow temperature range around 1 MeV. The nonlinear 
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corrections have a small effect on the corresponding shape of the "hump". The reason 
for this is a relatively high temperature (T ~ 1 MeV), at which the corresponding effects 
can occur. At such moderately high temperatures, the interplay between the two types 
of weak processes is substantially affected only if the nonlinearity (measured by SriQ/no) 
is well above 10%. 

In this study, we also found that the anharmonicity of density oscillations has 
an effect on the bulk viscosity, even though the effect was not large in the cases 
that we studied. For a strong anharmonicity, the bulk viscosity showed a substantial 
decrease. We also saw that different types of anharmonicity have slightly different 
qualitative as well as quantitative outcomes. This finding may suggest that some types 
of anharmonicity may be more efficient and, thus, lead to larger corrections to the bulk 
viscosity. In this preliminary study, we considered only a toy model to introduce the 
anharmonic oscillations. In the future, it may be interesting to study more realistic 
types of density oscillations that result from the actual nonlinear dynamics of stellar 
r-modes, produced by the gravitational emission. 

Here we studied the bulk viscosity only in the normal phase of dense quark matter. 
One may wonder, however, how the results will be modified if quark matter is a color 
superconductor [36] . As suggested by several existing studies in the linear regime, color 
superconductivity can have a large effect on the bulk viscosity [201 1211 EZl |38]. The 
additional effects due to the nonlinear regime may be much harder to predict. One of 
the complications comes from the fact that the weak rates in color superconducting 
matter have a strong dependence on the ratio of the superconducting energy gap 
and temperature, 4>/T. Moreover, in the suprathermal regime, their dependence on 
Sfii/T will most likely be nonpolynomial. These technical difficulties can be resolved in 
principle, but the outcome is not obvious and the corresponding study is outside the 
scope of the present paper. 
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Appendix A. Anharmonic oscillator 

Appendix A.l. Anharmonic oscillator with a cubic potential (Type I) 

The anharmonic oscillator with cubic potential U{x) = ttiUq + '^t) described by 
the following equation of motion: 



By making use of the known parametric solution to the above differential equation [39] 
and assuming that > is the maximum deviation from the equilibrium point x = 
in the positive x-direction, we find that the periodic solution exists for —1 < axo < 1/2. 




(A.l) 
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Figure Al. Solutions to the equation of motion of the anharmonic osciUator with 
a cubic potential for two values of the coupling constant. The inserts show the 
corresponding shapes of the potentials. 



(The apparent asymmetry between positive and negative values of ax^ is a result of the 
assumption that xq is the maximum deviation from the equilibrium point in the positive 
x-direction.) It is given in terms of the Jacobi elliptic function sn {u\m) as follows: 

'{t - to)uJo./a:^ ^ ^ 



x{t) = Xq 



3a_ 
2a*' 



-sn 



(A.2) 



where a± = ^ (l + 2a* ± + — 2a*)^ and a* = axo is a dimensionless 

parameter that measures the maximum deviation of the solution from the harmonic 
regime. Note that x(t) is periodic with the period given by 

4 



-K 



5. ^2 ^/ \3 

[a Xq) + —{a Xq) 



O [{axoy] ) , for axo^ 0,(A.3) 



2n ( 

where K {a_j a^^^ is the complete elliptic integral of the first kind. Two representative 
solutions are shown in Fig. lAll For anharmonic solutions of this type, the average 
kinetic energy is given by 



(-E'kin) = 



r, 



mx 



2^2 



dt 



TTIUqXq 



a Jo 



where the analytical expression for constant J-' reads 
4(1 + a* 



J" 



' a+ 



5a*a_ 




(A.4) 



(A.5) 



Here E (a_/a+) is the complete elliptic integral of the second kind. 
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Figure A2. Solutions to the equation of motion of the anharmonic osciUator with 
a quartic potential for two values of the coupling constant. The inserts show the 
corresponding shapes of the potentials. 



Appendix A. 2. Anharmonic oscillator with a quartic potential (Type II) 

The anharmonic oscillator with quartic potential U{x) = muQ + is described 
by the following equation of motion: 

X + co^x (1 + (3x^) = 0, (Type II). (A.6) 

Using the known parametric solution to the above differential equation [52] and assuming 
that Xo is the maximum deviation from the equilibrium point x = 0, we find that the 
periodic solution exists for /3* = /3xq > —1. The solution is given in terms of the Jacobi 
elliptic function, 



x{t) = XQ ^^t^^^?^* sd ( oooit - to) Vl+J^ 



This is a periodic solution with the period equal 

4 /3 



(3* 



2(1 + 



(A.7) 



271 



K 



2(1 + 13* 



UJQ 



1^4 + ^^P' 4 + {0^x1)]^ for /3 



0. 



Two representative solutions are shown in Fig. IA2[ For solutions of this type, the average 
kinetic can be given in the same form as in Eq. (lA.4p . but the value of the corresponding 
constant J-" is different. 



4(1 + 
3/3* 



E 



2(l+/3') 



K 



13* 



2(l+/3*) 



(A.9) 
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